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Timely overhauls produce in-flight service 
reliability greater than the calculated design 
reliabilities of the transmission components. 
Although necessary for aircraft safety, transmis- 
sion overhauls contribute to aircraft expense. 
Predictions of the transmission’s maintenance needs 
at the design stage should enable the development 
of more cost-effective and reliable transmissions 
in the future. 

This work estimates the frequency of overhaul 
and the number of transmissions or components 
needed to support the overhaul schedule. Two meth- 
ods based on the two-parameter Weibull statistical 
distribution for component life are used to esti- 
mate the time between overhauls. These methods 
predict transmission lives for maintenance sched- 
ules which (1) repair the transmission with a com- 
plete system replacement or (2) repair only failed 
components of the transmission . An example illus- 
trates the methods. 

INTRODUCTION 

Aircraft transmissions include bearings and 
gears which have finite fatigue lives with detect- 
able failure warnings. The two-parameter Weibull 
distribution statistically describes the drive 
system bearing and gear life [1-4]. 

The in-flight service re 1 i ability of aircraft 
transmissions is much higher than the design reli- 
ability of their components. Transmission over- 
hauls provide the difference. By monitoring the 
onset of potential fatigue failures, one can use 
just-in-time overhauls to maintain the transmission 
economically and reliably [5]. 

A two-parameter Weibull distribution provides 
the transmission system life model for repairing a 
transmission with full-system replacement [6,7]. 

The sum of the component failure rates predicts the 
repair frequency for maintenance with partial- 
replacement repair [ 8,9 ] . 

Renewal theory Is a secondary statistical 
model that describes the maintenance process. It 
estimates the number of replacements needed to 
maintain transmission reliability with a specified 
maintenance schedule. The theory considers the 
ongoing sequence of use, failure onset, repair, and 
return to use. For this sequence, renewal theory 
predicts the frequency of component replacement and 
the number of replacements needed to support the 
service maintenance schedule [ 10-12 ] , 


Confidence theory complements these statistics 
with estimates of the likelihood of the predic- 
tions. Higher confidence levels require more 
spare parts to cover a greater range of possible 
situations [ 11 , 12 ] . 

This work presents these theories and applies 
them to a simple transmission (Fig. 1) to show 
their use. Estimates of drive system component 
failure onset rate and replacement needs are essen- 
tial in design. They allow one to compare the 
worth of different designs from a maintenance cost 
perspective, and they help assess the cost of oper- 
ating a proposed drive system design. 

COMPONENT RELIABILITY 

The two-parameter Weibull distribution is a 
statistical function commonly used to describe 
fatigue life data. It can describe a variety of 
life patterns in which the reliability of a compo- 
nent is the complement of its probability of 
failure. 



Fig. 1. Single-mesh transmission. 




In statistics, reliability is a double nega- 
tive. Reliability, or the act of surviving, is the 
state of not having failed. Statistics count sin- 
gle events such as the act of failing. A part can 
fail only once, whereas it survives for its entire 
life. The probability of failure for the two- 
parameter Weibull distribution, which is a direct 
statistic, is 

F - 1.0 - e' (,/e)b - 1.0 - R 0) 


either equation generates a straight-line plot as 
shown in Fig. 2. The plot is a probability graph 
for the two-parameter Weibull distribution. 

The average life is the mean time to failure 
(MTTF), which is the sum of all times to failure 
divided by the total number of failures. For a 
continuous distribution, the total number of fail- 
ures is unity, and the sum of all lives to failure 
is the integral of time or life times the probabil- 
ity density function. Integrating from zero to 
infinity gives the mean life: 


The derivative of Eq. (1) with respect to life 
is the probability density function f: 


f 


b-l 


-(*/<?)> 


( 2 ) 


The probability density function is a histo- 
gram of life failures for a unit population. The 
Weibull reliability function can be expressed as a 
log reciprocal: 
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In working with the high-reliability range, 
? 10 1 ife often replaces the characteristic 
8 as the scaling parameter. In terms of 
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Even though it is cumbersome, manufacturers 
use Eq. (4) as the two-parameter Weibull distribu- 
tion of bearings to place 90-percent reliability 
lives in the catalogs [13!] . 

In both Eqs. (3) and (4), the log of the reli- 
ability reciprocal is proportional to the life 
raised to the Weibull slope. Taking the log of 
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Fig. 2. Two-parameter Weibull probability plot. 
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Substituting the probability density function 
of Eq. (2) for the two-parameter Weibull distribu- 
tion and integrating yields the well-known gamma 
function T multiplied by the characteristic 
life 8: 


l . MTTF = 8 f 


The solid curve in Fig. 3 is a plot of the 
ratio of the two-parameter Weibull mean life to the 
characteristic life versus Weibull slope. The mean 
life equals the characteristic life at b = 1.0, 
drops below the character! st ic life to a minimum 
relative value at b = 2.15, and then increases 
back to the characteristic life as h approaches 
infinity. When b is infinite, the distribution 
is an impulse with all lives equal to the charac- 
teristic life. 

By a similar integration, one can find the 
standard deviation of the two-parameter Weibull 
distribution. The standard deviation is the square 
root of the second moment of the component life 
distribution about the mean. 
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Fig. 3. Average life and standard deviation of life 
ratios to characteristic life for a Weibull distribution 
as a function of the Weibull slope. 
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In terms of the gamma function, the standard devia- 
tion of the two-parameter Weibull distribution is 




( 12 ) 



The standard deviation of a distribution is a 
measure of the scatter of the distribution. It is 
valuable in estimating a confidence limit for the 
average life. 

The broken curve in Fig. 3 is a plot of 
the ratio of the standard deviation of the two- 
parameter Weibull distribution to its characteris- 
tic life versus the Weibull slope. At a slope of 
one, the distribution is the exponential distri- 
bution, which has a large scatter. As the slope 
increases to two, the scatter decreases rapidly and 
continues to decrease with increasing slope. 

SYSTEM LIFE WITH FULL REPLACEMENT 

To model the transmission life based on full 
replacement, one must have a life model for the 
transmission as a complete system which treats the 
system as a single component. The life of a drive 
system can be considered to be a strict series 
probability model of the lives of its components 
[7]. In this model, the reliability of the system 
R s is the product of the reliabilities of all the 
components : 

r s - n Ri 


The high speed of drive system components and 
the spray of loose debris warrant the strict series 
probability model. If any component fails, debris 
may be present which could accelerate the fatigue 
damage in other components. Therefore, the drive 
system will need an overhaul to return it to a high 
state of reliability when any element fails. 

The log of the reciprocal of Eq. (9) is 



The straight-line reliability relationship of 
Eq. (12) can be fit numerically to the more exact 
relationship of Eq. (11) with a linear regression. 
The slope of the fitted straight line is the drive 
system Weibull slope b $t and ? 1Q s is the life at 

which the drive system reliability R $ equals 
90 percent on the straight line. 

SYSTEM LIFE WITH PARTIAL REPLACEMENT 

To model the transmission live based on par- 
tial replacement, one can treat the full system as 
a collection of independent components. Separate 
analysis of each component will predict the number 
of replacements needed. If no two components are 
repaired at the same overhaul, the maximum number 
of overhauls is equal to the sum of all individual 
component replacements. One can estimate the mean 
time between overhauls as the total service time 
divided by the number of replacements for this 
component sum repair calculation and for the full- 
system repair. 

RENEWAL THEORY 

Renewal theory estimates the number of 
replacements as a function of the component failure 
distribution and its life [ 10-12 ] . It assumes that 
failed components will be replaced just before they 
fail, which models an unending sequence of use and 
repair. Aircraft drive system maintenance follows 
this pattern closely. 

The mean number of failures is the infinite 
sum of the probabi 1 i ties of at least i failures 
in the life period t. This function, M(*), is the 
renewal function. It is expressed as 

M(0 = Fj(i) * M(< - x)f(x) dx (13) 

The derivative of the renewal function with 
respect to life is the renewal density function: 



m(t) - fi(«) m(f -x)f(x) dx (14) 


and substitution of Eq. (4) into Eq. (10) for each 
component yields 




(id 


In Eq. (11), ? s is the life of the entire drive 
system for the system reliability R s . It is also 
the life of each component at the same drive system 
reliability R . For consistency in Eq. (11), all 
the component lives must have the same counting 
base of hours. 

Equation (11) is a two-parameter Weibull dis- 
tribution only when all the Weibull exponents b i 
are equal. However, a two-parameter Weibull dis- 
tribution can approximate Eq. (11) quite well. 

Equation (12) is the drive system two- 
parameter Weibull relationship. It includes the 
system reliability parameters b $ and t 10 s 4 


These equations give the number of replace- 
ments needed to support a maintenance schedule. 
Their solution involves a series of convolution 
integrals that can be performed on any failure 
distribution. However, the solution, which is an 
oscillation of replacement numbers about a straight 
line, is not easily obtained. The solid curve of 
Fig. 4 shows the renewal function for a component 
with a two-parameter Weibull reliability, 

0 - 5000 hr and b = 1.5. Tabulated solutions to 
the renewal function for the two-parameter Weibull 
distribution are available [12]. 

An approximation for the renewal function [11] 
is 

2 2 

M e (t ) . J_ -Ifi _ gf (15) 
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The accuracy of this approximation increases 
as « increases. Equation (15) is an asymptote 
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for the true renewal function of low-scatter dis- 
tributions. For high-scatter di stributions , it 
approximates the true renewal function closely. 

The renewal function is the probability of 
replacement for a single component. Its value goes 
above one because multiple replacements can occur. 
For a set of Q identical components, the total 
number of replacements is the product: 

N r -QM(|) (16) 


Estimates of replacement inventory need a 
margin for variations from the mean, as do repair 
frequency estimates. Confidence statistics based 
on the renewal standard deviation provide one means 
for determining this margin. The broken curve of 
Fig. 4 is a plot of the renewal function standard 
deviation versus life for the component with a 
characteristic life of 5000 hr and a Weibull slope 
of 1,5 for which the solid curve of Fig. 4 plots 
the renewal mean. 

The approximation for the standard deviation 
of the renewal function uses the third moment of 
the life distribution . For the two-parameter 
Weibull distribution, the third moment is 


The standard deviation of the renewal function 
gives a measure of the scatter in replacement needs 
from one sample to the next. Estimates of replace- 
ment inventory need a margin for variations from 
the mean, as repair frequency estimates do. Confi- 
dence statistics provide one means for determining 
this margin. 

CONFIDENCE STATISTICS 

In predicting replacement rates and mainte- 
nance inventories, direct theory provides mean or 
• average’ estimates. These estimates come from the 
statistics of a universal population. With enough 
cases, they will be the true average values. 

In any real situation, the number of drive 
systems under service is a limited sample. Confi- 
dence statistics estimate how differently a small 
sample may behave from its universal population. 

It uses the standard deviation of the universal 
failure distribution and the sample size to esti- 
mate the mean of the sample. Confidence intervals 
are shown by the broken lines in Fig. 2, 

For many samples of the same size, the mean of 
the samples has a normal distribution about the 
overall mean. The standard deviation of the means 
is 
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Figure 5 shows the third moment of the two- 
parameter Weibull distribution divided by the cube 
of the characteristic life versus the Weibull 
slope. 

The approximat ion for the standard deviation 
of the renewal function is [U] 
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where Q is the size of the sample. 

In reliability predictions, the lower confi- 
dence bound is valuable in aircraft applications. 
This confidence distribution estimates the life at 
which a large percentage of the samples of a given 
set will survive. For a high confidence, this life 
is less than the mean life for the entire popula- 
tion. For a 90-percent confidence, 

*2v ( 90 = f av “ z 10 cr av (20) 




Fig. 4. Renewal function and renewal function 
standard deviation for a two-parameter Weibull 
distribution with 0 = 5000 hr, and B=1.5. 
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Fig. 5. Third life moment to characteristic life cubed 
ratio for a Weibull distribution as a function of the 
Weibull slope. 


Table 1. - Single-mesh transmission properties 



Load, 

kN 

10-percent 

life, 

6*10' 

10 b cycles 

10-percent 

life, 

* 10' 
hr 

Average 

life, 

* av ' 

hr 

Standard 
deviation 
of life, 

hr 

Bearing 1 

4.1 

317 

2640 

16 187 

13 570 

Bearing 2 

10.3 

578 

4820 

29 554 

24 776 

Pinion 

5.8 

298 

2486 

5 426 

2 324 

Bearing 3 

10.3 

868 

7230 

44 330 

24 460 

Bearing 4 

4.1 

475 

3960 

24 280 

20 097 

Gear 

5.8 

380 

3170 

6 920 

2 963 

Transmi ssion 


127 

1060 

3 990 

2 600 


Since the behavior of samples differs from the 
behavior of the "ideal* distribution, confidence 
estimates help one to see the effects of sample 
size on the life and replacement estimates. 

EXAMPLE 

For the single-mesh transmission shown in 
Fig. 1, the 90-percent reliability lives for the 
bearings and gears are shown in Table 1. The 
Wei bull slope for the bearings is 1.2, and for the 
gears is 2.5. For a fleet of Q = 50 aircraft, we 
would like to estimate the number of overhauls in 
the first 10 000 hr of service and the number of 
replacement components needed to support these 
overhauls. Two types of overhaul are treated - 
full replacement and f ai led-component replacement 
only. All estimates will be with 90-percent confi- 
dence for the 50 aircraft sample size. From 
Eq. (6), the average lives were determined for each 
component. From Eq. (8), the standard deviations 
for each component were determined. The results 
are shown in Table 1 . 

Full Replacement 

To treat the transmission as a complete system 
undergoing ful 1 -repl acement repair, one can use the 
two-parameter Weibull system model of Eq. (12). 

The parameters b $ and t lQ s for Eq. (12) come 
from a least squares fit to'Eq. (11). The two- 
parameter Weibull slope is b s = 1.57 for the 
transmission, and the system 90-percent reliability 
life is * 10s = 1060 hr. From Eq. (3), the trans- 
mission characteristic life is 6 = 4440 hr. For 
these data, the transmission average life is 
f av - 3990 hr with a standard deviation of 
a f = 2600 hr. 

The renewal function can estimate the number 
of transmissions needed for full replacement in a 
continual sequence of failure warning, repair, and 
return to service for the 50 aircraft. For an 
average life of 3990 hr and a standard deviation 
life of 2600 hr, Eqs. (15) and (16) give the total 
number of replacements in the period from 0 to 
i hr. 

From Eq. (17), the third moment of the trans- 
mission life distribution is /< 3 = 1.62xlO u hr 3 , 
and Eq. (18) gives the standard deviation of the 
renewal function for the transmission. 

Equations (18) and (19) give the standard 
deviation of the number of replacements in the 
period from 0 to {. Finally, a relationship simi- 
lar to Eq. (20) gives the replacement estimate for 
complete transmissions with a 90-percent confidence 
that the replacements will be less. 


For the first 10 000 hr of operation, this 
procedure estimates an average number of 111 
replacements for the 50 aircraft. A confidence 
limit of 90 percent boosts this estimate to 121 
transmission replacements for 500 000 fleet service 
hours. This represents a mean time between over- 
hauls of 4130 hr and a total spare parts require- 
ment of 726 parts. 

Partial Replacement 

When only the f ai 1 i ng components are replaced, 
the renewal function can estimate the number of 
replacements needed, also. Applying the calcula- 
tions of this procedure for each of the six compo- 
nents in the transmission estimates the number of 
components needed to support a partial -repai r main- 
tenance schedule with a 90-percent confidence. 

Table 2 summarizes these calculations for the 
four bearings and two gears in the transmission. 
Adding the total number of components that renewal 
theory estimates will need to be replaced yields 72 
bearings and 135 gears, for a total of 207 compo- 
nents. This total of 207 spare parts is signifi- 
cantly less than the 726 parts required by the 121 
full -transmi ssion replacements required of the 
other service procedure. 

If each component failure required its own 
overhaul, then 207 overhauls would be required with 
the same 90-percent confidence as used for the 
ful 1 -repl acement calculations. Dividing the 
500 000 fleet service hours by the maximum number 
of 207 overhauls yields an estimate for the mean 
time between overhauls equal to 2420 hr. This is 
1710 hr less than the mean time between overhauls 
for ful 1 -transmi ssion replacement because it does 
not consider repair of components near failure in a 
maintenance session. 

By only replacing the failed components, one 
would need 86 more overhauls, but 519 less parts. 


Table 2. - 10 000-hr repair estimates 
for 50 transmissions 


Transmission 

component 

Replacements 

requi red 

Mean 

Standard 

deviation 

90-percent 

confidence 

Bearing 1 

23 

5.2 

30 

Bearing 2 

9 

4.1 

15 

Pinion 

72 

4.7 

78 

Bearing 3 

4 

3.6 

8 

Bearing 4 

13 

4.4 

19 

Gear 

52 

4.3 

57 

Total 

173 


207 
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The same high reliability would be present for both 
maintenance procedures because of the on-board 
failure monitoring system. These estimates are for 
cost and scheduling purposes only. 

SUMMARY OF RESULTS 

Two methods of estimating the time between 
transmission overhauls and the number of replace- 
ment components needed are presented. The first 
treats full replacement of failed transmissions, 
whereas the second treats replacement of failed 
components only. Confidence statistics are applied 
to both methods to improve the statistical estimate 
of sample behavior. 

The method to predict overhaul timing with 
full replacement is based on a two-parameter 
Wei bull system life model. The relationship 
between the system life model and the component 
life models is presented. In addition, formulas 
for the mean and standard deviation of the two- 
parameter Wei bul 1 distribution are given. 

Renewal theory is presented as a tool to esti- 
mate the number of replacements in a transmission 
undergoing a consistent maintenance procedure. The 
theory is useful for estimating replacements for 
both full and partial transmission-replacement 
procedures. Approximation formulas are given for 
the mean and standard deviations of the renewal 
function. These approximations are valid for two- 
parameter Weibull distribution lives amongst 
others. Formulas for sample replacement numbers 
are given in terms of the renewal function. 

Single-sided confidence theory is presented 
forthe replacement number and overhaul timing 
estimates. A transmission example is presented to 
illustrate the methods. Comparisons of overhaul 
timing and spare-part requirements are made in the 
example between ful 1 -transmi ssion replacement and 
partial component-replacement overhauls. High 
reliability is assured for the transmi ssions in 
both cases by the on-board monitoring system. 

NOMENCLATURE 
b Weibull slope 

e base of the natural log 

F probability distribution function (probabil- 

ity of failure) 

f probability density function 

Ln natural log 

? 1 i f e , h r 

M renewal function 

M e approximate renewal function 

MTTF mean time to failure 

m renewal density function 

N r number of replacements 

Q sample size 

R reliability (probability of survival) 

x integration time variable, hr 

z 10 number of standard deviations from the mean 

which cuts off a IG-percent population tail 

T gamma function 

6 characteristic life, hr 

/< 3 third moment of a probability density 

function 


a standard deviation 

o f standard deviation of Weibull function 

a me standard deviation of renewal function 

Subscri pts : 
av average or mean 

i index 

n number of components 

s system 

1 index value 

10 10-percent failure, 90-percent 

reliability 

90 90-percent confidence 
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